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Abstract 

Although the study of weak convergence of superposition of point processes to the 
Poisson process dates back to the work of Grigelionis in 1963, it was only recently that 
Schuhmacher (2005a) obtained error bounds for the weak convergence. Schuhmacher 
considered dependent superposition, truncated the individual point processes to 0- 
1 point processes and then applied Stein's method to the latter. In this paper we 
take a different approach to the problem by using Palm theory and Stein's method, 
thereby expressing the error bounds in terms of the mean measures of the individual 
point processes, which is not possible by Schuhmacher's approach. We consider locally 
dependent superposition as a generalization of the locally dependent point process 
introduced in Chen and Xia (2004) and apply the main theorem to the superposition 
of thinned point processes and of renewal processes. 

Key words and phrases. Stein's method, Poisson process approximation, dependent su- 
perposition of point processes, sparse point processes, thinned point processes, renewal 
processes. 
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1 Introduction 

The study of weak convergence of superposition of point processes dates back to Grigelionis 
(1963) who proved that the superposition of independent sparse point processes converges 
weakly to a Poisson process on the carrier space IR + . His result was subsequently extended 
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to more general carrier spaces by Goldman (1967) and Jagers (1972) [see Cinlar (1973) and 
Brown (1978) for further discussions]. It was further extended to superposition of depen- 
dent sparse point processes by Banis (1975, 1985), Kallenberg (1975), Brown (1979) and 
Banys (1980). For a systematic account of these developments, see Kallenberg (1983). 

Surprisingly, it was only recently that error bounds for such convergence of point pro- 
cesses were studied. Using Stein's method for Poisson process approximation as developed 
by Barbour (1988) and Barbour and Brown (1992), Schuhmacher (2005a) obtained an er- 
ror bound on the d,2 Wasserstein distance between a sum of weakly dependent sparse point 
processes {£ n j, 1 < i < /c n }neN and an approximating Poisson process. As he truncated the 
sparse point processes to 0-1 point processes as in the proof of the Grigelionis theorem, his 
error bound contains the term YliZi ^[£ni(B) > 2], whose convergence to for every bounded 
Borel subset B of the carrier space is a condition for the Grigelionis theorem to hold. A con- 
sequence of such truncation is that the mean measure of the approximating Poisson process 
is not equal to the sum of the mean measures of the individual point processes. 

In this paper we take a different approach to the Poisson process approximation in which 
we do not use the truncation, but apply the Palm theory and express the error bounds in 
terms of the mean measures of the individual sparse point processes. Such an approach also 
ensures that the mean measure of the approximating Poisson process is equal to the sum of 
the mean measures of the sparse point processes. 

As in Schuhmacher (2005a), we study the dependent superposition of sparse point pro- 
cesses. But we consider only locally dependent superposition which is a natural extension 
of the point processes ^IiS^ studied in Section 4 of Chen and Xia (2004), where S x is the 
point 's are 5-valued independent random elements with S a locally compact 

metric space, the indicators Zj's are locally dependent and the Jj's are independent of the 
tfi's. 

In our main theorem (Theorem 2.1), with the help of Lemma 3.1 in Brown, Weinberg 
and Xia (2000), it is possible to recover a factor of order 1/A from the term 1/(|E]W| + 1). 
Hence the error bound on the d,2 Wasserstein distance yields the so-called Stein factor 1/A, 
by which approximation remains good for large A, a feature always sought after for Poisson- 
type approximations. In the error bound obtained by Schuhmacher (2005a), a leading term 
does not have the Stein factor [see Remark 4.4 for further details]. 

Our main theorem and some corollaries are presented in Section 2. Applications to 
thinned point processes and renewal processes are given in Sections 3 and 4 respectively. 

2 The main theorem 

Throughout this paper, we assume that T is a locally compact metric space with metric do 
bounded by 1. In estimating the error of Poisson process approximation to the superposition 
of dependent point processes {Sj, i 6 1} on the carrier space T with X a finite or countably 
infinite index set, one natural approach is to partition the index set X into {{i}, Xf, X™}, 
where If is the set of the indices of the point processes which are strongly dependent of 
Sj and Tf the set of the indices of the point processes which are weakly dependent of 
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[see Schuhmacher (2005a)]. Another approach is to divide the index set according to various 
levels of local dependence, a successful structure for studying normal approximation [see 
Chen and Shao (2004)]. The later approach has been generalized by Barbour and Xia (2006) 
to randomly indexed sums with a particular interest in random variables resulting from 
integrating a random field with respect to a point process. 

Parallel to the local dependence structures defined in Chen and Shao (2004), we intro- 
duce: 



[LD1] For each % G X, there exists a neighbourhood Ai such that i G Ai and Sj is independent 
of {5,, j G A^}. 

[LD2] Condition [LD1] and for each i G X, there exists a neighbourhood Bi such that 
Ai C Bi and {Ej, j G A { } is independent of {Ej, i G Bf\. 



The index set X in [LD1] and [LD2] will be assumed to be finite or countably infinite in 
this paper, although it may be as general as that considered in Barbour and Xia (2006). The 
superposition of {Sj : % G X} which satisfies the condition [LD1] is more general than point 
processes of the form Yl^^Ui where the I^s are locally dependent indicators with one level 
of dependent neighborhoods in X (that is, the lis satisfy [LD1] in Chen and Shao (2004, 
p. 1986)). Such a point process is a typical example of locally dependent point processes 
defined in Chen and Xia (2004, p. 2548). Likewise, the superposition of {Ej : i G X} which 
satisfies the condition [LD2] is more general than point processes of the form ^2 IiSu i where 
the Ij's are locally dependent indicators with two levels of dependent neighborhoods in X 
(that is, the I^s satisfy [LD2] in Chen and Shao (2004, p. 1986)). 

Three metrics will be used to describe the accuracy of Poisson process approximation: 
the total variation metric for Poisson random variable approximation d tv , the total variation 
metric for Poisson process approximation dxv an d a Wasserstein metric d 2 [see Barbour, 
Hoist and Janson (1992) or Xia (2005)]. 

To briefly define these metrics, let "H be the space of all finite point process configurations 
on r, i.e., each £ G "H is a non-negative integer-valued finite measure on T. Let K. stand for 
the set of c?o-Lipschitz functions k : r — > [—1, 1] such that | k(a) — k((3) \< do(a,(3) for all 
a, (3 G T. The first Wasserstein metric d\ on "H is defined by 

r o, if i6i = i6i = o, 

d 1 (£ 1 ,£ 2 ) = \ 1, if |6| ^ |6|, 

( \Zi\- 1 sup kelc \fkdti 1 -fkd{i 2 \, if |6| = |6I >o, 

where |6| is the total mass of 6- A metric d[ equivalent to d\ can be defined as follows [see 
Brown and Xia (1995)]: for two configurations 6 = Y^i=i anc ^ £2 = YlT=i with m > n, 



n 



G?i(6, 6) = min N d (yi, z n(i) ) + (m 



i=l 



n) 



where tt ranges over all permutations of (1, . . . , m). Both d\ and d[ generate the weak topol- 
ogy on % [see Xia (2005, Proposition 4.2)] and we use B(H) to stand for the Borel a- 
algebra generated by the weak topology. Define three subsets of real valued functions on % 
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as T tv = {l A m : A C Z+}, JF dl = {f : |/(&) - /(&)| < for all 6, 6 G 

and J-tv = : ^4 G B(H)}. Then the pseudometric <i to and the metrics d 2 and cItv are 

defined on probability measures on "H by 



^(Qi, Q 2 ) 

4(Qi,q 2 : 

<W(Ql, Q 2 ) 



inf PflXxl ^ |X 2 |) = sup 
inf EfdifXi,^)] = sup 
inf P(Xi ^ X 2 ) = sup 



/dQi- 



where the infima are taken over all couplings of (Xi, X2) such that £(Xj) = Qj, z = 1, 2 and 
the second equations are due to the duality theorem [see Rachev (1991, p. 168)]. 

To bound the error of Poisson process approximation, we need the Palm distributions 
Q a of a point process X 2 with respect to a point process X\ with finite mean measure v at 
a. When X\ is a simple point process, i.e., it has at most one point at each location, the 
Palm distribution Q a may be intuitively interpreted as the conditional distribution of X2 
given that X\ has a point at a. More precisely, let B{T) denote the Borel a— algebra in V 
generated by the metric d and define the Campbell measure C of (Xi, X 2 ) on B(T) x B(H): 

C(B x M) = EfXx^lj^M], B G B{T), M G B(H). 

Since the mean measure v oiX\ is finite, by [Kallenberg (1983, 15.3.3)], there exist probability 
measures Q a on B(H) such that 



E[Xi(S)lx ae jif] 



Q a (M)u(da), V 5 G S(r), M G 



(2.1) 



which is equivalent to 



Qa{M) 



'E{l [ x 2e M]X 1 (da)} 
v(da) 



Vie B{H), aeT i/-a.s. 



[see Kallenberg (1983, section 10.1)]. It is possible to realize a family of point processes Y a 
on some probability space such that Y Q ~ Qa, and we say that Y a is a Palm process of X 2 
with respect to Xi at a. Moreover, when X\ = X 2 , we call the point process Y a — 5 a the 
reduced Palm process of X 2 at a [see Kallenberg (1983, Lemma 10.2)]. 

As noted in Goldstein and Xia (2006), when T is reduced to one point only, the Palm 
distribution of X 2 (with respect to itself) is the same as the size-biased distribution and 
general guidelines for the construction of size-biased variables are investigated in Goldstein 
and Rinott (1996). 



Theorem 2.1 Let {Hj, j 6 1} be a collection of point processes on T with mean measures 
Xi, i G X respectively. Set H = w ^ ^ean measure denoted by X and assume that 
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A := A(r) < oo. If [LD1] holds, then 

4,(£(S),Po(A)) < ^T^EV f {\\m-\V i>a \\ + 11^1-1-^11} Xi(da), (2.2) 

A id J r 

d 2 (£(~),Po(A)) < E^^ + p^) jfdiC^V^Cda) 

+ E (t + E |SW| + i ) ^ J T d i(^ E ^)Xt(da), (2.3) 
<W(£(H),Po(A)) < EV /{||^-^ )a || + 11^-5^11} \(da), (2.4) 

where = ^ 6j4l? Sj ; V, = 5^ J - eJ 4 j \{j} ^ e reduced Palm process o/Sj at a, Vj )Qi 

£ne Pa/m process of Vi with respect to S, at a suc/i £/iat -f Vi )Q + Sj )Q + 5 a tne 
Pa/m process of S respect to Sj at a, ana 1 || • || denotes the variation norm of signed 
measure. Under the condition of [LD2], (2.2) and (2.4) remain the same but (2.3) can be 
further reduced to 

d 2 (£(H),Po(A)) < W^ + E= |E / d[(V t ,V ha )X t (da) 



(2.5) 



< 



where A.j = Aj(r) ana 1 



+E(|^| • 15,1) + A 4 2 + E (|Hi| 2 ) - A,}, (2.6) 



Proof. We employ Stein's method for Poisson process approximation established in Bar- 
bour (1988) and Barbour and Brown (1992) to prove the theorem. To this end, for suitable 
measurable function h on H, let 

= Jm + S a ) - h(£)]\(da) + j[h{i - 8 X ) - h(OWx). 

Then A defines a generator of the spatial immigration-death process with immigration in- 
tensity A and with unit per capita death rate, and the equilibrium distribution of the spatial 
immigration-death process is Po(A) [see Xia (2005, section 3.2) for more details]. The Stein 
equation based on A is 

Ah(t) = f(t)-Po(\)(f) (2.7) 



5 



with solution 



hf(0 



[E/(Z c (t))-Po(A)(/)]dt, 



where {Z^{t), t > 0} is the spatial immigration-death process with generator A and initial 
configuration Zg(0) = £. To obtain bounds on the errors in the approximation, we need to 
define 



Ah f (Z;x) 
A 2 h f (£;x,y) 
A 2 h f (Z,r);x) 



hf(t + 6 x )-h f (t), 
Ah f (Z + 6 x ;y)-Ah f (t;y), 
Ah f (£;x) - Ah f (ri;x), 



for corresponding test functions /. Xia (2005, Propositions 5.6 and 5.12) [see also Barbour 
and Eagleson (1983), Barbour and Brown (1992)] and Xia (2005, Lemma 5.26) state that, 
for all x, y G T, 



\A 2 h f (^x,y)\<^-^, V/eJ to ; 
\A 2 h f (C;x,y)\<l, V/eJ rv ; 



(2.8) 
(2.9) 



|A 2 M^;x)| < (x + \ v \ A 2 |jj + l ) d '^), (2.10) 



Now, since + y i a + E iyOL + 5 a is the Palm process of 5 with respect to Hj at a, it 
follows from (2.1) that 



E J [h(E) - h(E - 6 a )]E{da) 

Ve/ [h(E)-h(E-6 a )]Ei{da) 
iex Jr 

Ve/ Ah(E^ + V ha + E iia ;a)Xi(da). 



On the other hand, by the Stein equation (2.7), we have 
|E/(H)-Pb(A)(/)| 



= mj[hf(E + 8 a ) - h f (E)]X(da) + Bj[h f (E - 8 X ) - h f (E)]E(dx) 

= E / i Ah f& a) - Ahf(E® + V i>a + E i>a ; a)}Xi(dc 

< E / i\ Ah f(Z {i) +Vi + E t ; a) - Ah f (E® + V i>a + E,; a)\ 

Ah f (E® + V ha + E i>a - a)|}A 4 (da) 



\Ah f {E^ + V ha + E l]a ) 



(2.11) 



(2.12) 
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To prove (2.2), we note that the test functions / G Ttv satisfy /(£) = f(\C,\S z ) for a fixed 
point z G T, and so we have = Hence, for all 77, 6, 6 G 

|A/ i/ (77 + a;a)-A/i / (7 7 + 6;a)| 

= |AM»7 + (161 V 161)5,; a) - Ah f ( V + (|6| A |6|)^; a)\ 
II6H6H 

< ^ lA^ + deilA^I+i-l)^;^^! 
i=i 

< 1161- 1611 (2.13) 

where the last inequality is due to (2.8). Combining (2.13) with (2.12) yields (2.2). 
Next, (2.10) and (2.11) imply that for / G 

|E/(H)-Pb(A)(/)| 

Because d'^Vi+Ei, Vi )0l +Ei j0l ) < d[(Vi, Vi^+d'^Ei, E ij(X ) and, for each i G X, S; is independent 
of E)W, (2.3) follows. On the other hand, due to the independence between {Vi,Sj} and 
{Hj, j G 5f} implied by [LD2], (2.5) is immediate. To prove (2.6), one can verify that 

Var(X;|H,-|) = £ £ covdS^U^I) 
yes? / jieflfiaeBfnA^ 

and that E YljeB? I^j'I = SjeB c "\r Hence (2.6) follows from Lemma 3.1 in Brown, Weinberg 
and Xia (2000) and the fact that d'^V^V^a) < \Vi\ + \V it0l \ and d[(Ei, Hj )Ct ) < |Sj| + |H i)Ct |. 

Finally, we show (2.4). For £1, £2 G we define 



£1 A £ 2 = J^(ay Aa 2j )4 i , 



where {xi, . . . ,Xk} is the support of the point measure 6 + 6 so that 6 = Y^j=i a ij&xj for 
z = 1,2 with o^'s being nonnegative integers. Then, for all / G Ttvi V, 6> £2 G "H, 

^(77 + 6; a) - A/i/fa + foaOl 

< |A/i / (77 + 6;a) - A/i / (77 + Ci A£ 2 ;a)| + |A/i/(t7 + 6; a) - A/i/fo + 6 A 6; a)| 

< (161 -16 a 61) + (161 -16 a 61) 

= 116-611, (2.14) 
where the last inequality is due to (2.9). Applying (2.14) in (2.12), we obtain (2.4). ■ 
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Corollary 2.2 With the notation in Theorem 2.1, i/{Sj, i £ X} are a// independent, then 

1 - e" A 



4,(£(s),Po(A)) < 



A 



- E E 



iex 



|S iiQ || Xi(da) 



d 2 (£(E),Po(A)) < + E ^ )E / <(5 i ,5 i , a )A i (da 

^/k(1 + k/4) + 1 + k/2 



1^ + 2.5- 



iex 



<W(£(3),Po(A)) < EV / ll^-H^llA^a 



(2.15) 
(2.16) 



(2.17) 
(2.18) 



where k 



A— maxj^x ^j+l * 



Proof. Let A, = B { = {i}, then (2.15), (2.16), (2.17) and (2.18) follow from (2.2), (2.5), 
(2.6) and (2.4) respectively. ■ 

Corollary 2.3 [cf Chen and Xia (2004, Theorem 4.1)] Let {ij, i £ X} 6e dependent indica- 
tors withX a finite or countably infinite index set and let {Ui, i £ X} be V -valued independent 
random elements independent of {/j, z £ X}. Define H = X^ex-^4 mean measure X, 
let E/j = and assume that A = ^2 ieI Pi < 00 ■ -^ 9r eac ^ * ^ ^ ^« ^ e ^ e se ^ °/ indices 
of those Ij 's which are dependent on Ii, that is, Ii is independent of {Ij : j £ A?}. Tnen 



4,(£(S),Po(A)) < 



1-e 
A~~ 



Ei E w, + E^- 



iex lieAA{j} 



d a (£(3),Po(A)) < E£ E (t + ^TI 

iex jeAA{i} V 



iex je^i 



3.5 



££(t+ e 



2.5 



S'i + l 



PiPj, 



d TV (c(E),Po(x)) < El E mij+^PiPj 

iex I jeAA{i} ieAi 



(2.19) 



(2.20) 



(2.21) 



where Si = Y^jgA-^j- F° r ea °h i £ X, fei 6e £ne se£ of indices of those Ii's which are 
dependent on {Ij, j £ A{\ so that {Ij : j £ A{\ is independent of {// : I £ -B?}. Tnen 
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(2.19) and (2.21) remain the same but (2.20) can be further reduced to 
d 2 (£(H),Po(A)) < E E ( 3 4 + ^WTl) m ^ ) 



< 



*ez jeAi\{i} 

('■ ._ 

iei V A A/j 6 flfPj + 1 / \iGMW 

(2.23) 

where Wi = J2j^Bi h an( ^ 



L J2, 



Proof. Set Sj = Ij^, i E Z, then is independent of {Hj : j ^ Aj}, so [LD1] holds, E ijCe = 0, 
and the claims (2.19), (2.20) and (2.21) follow from (2.2), (2.3) and (2.4) respectively. On 
the other hand, {Ej : j E Ai} is independent of {Ej : j E" Bi}, so [LD2] holds and (2.22) 
and (2.23) are direct consequences of (2.5) and (2.6). ■ 

A typical example of Poisson process approximation is that of the Bernoulli process 
defined as follows [see Xia (2005, section 6.1) for further discussion]. Let I\, ■ ■ ■ , I n be 
independent indicators with 

P(/ t = 1) = 1 - = 0) = Pi , % = 1, ■ • • ,n. 

Let r = [0, 1], E — Yli=i^i/n an d A = Yli=iPA/n be the mean measure of E. We set 
Hj = IiSi/ n , i — 1, • • • , n, then the reduced Palm process of at a E V is E it0l = and the 
Palm distribution of Ej with respect to point process X; at a for j ^ i is the same that of 
Ej. Hence, Corollary 2.2 together with (2.16) and Chen and Xia (2004, Proposition 4.5) can 
be used to obtain immediately the following (known) result. 

Example 2.4 [Xia (2005, section 6.1)] For the Bernoulli process E on T = [0, 1] with mean 
measure A, 

l-e- A ^ 2 



cZ to (£(H),Po(A)) < J^pI 



i=l 



d TV (£(E),Po(X)) < £ 



i=l 

fi — 

2 



d 2 (£(H),Po(A)) < Yp: 

— i=i 
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Example 2.5 Throw n points uniformly and independently onto the interval [0, n] and let 
S be the configuration of the points on [0, T] := T with n ^> T and A be the mean measure 
of S, then 

d 2 (£(£),Po(A))<-^-. 



Proof. Let 7j = 1 if the zth point is in T and if it is not in T. Then the configuration of 
the ith point on T can be written as Sj = and H = Ym=i where Wj's are independent 
and identically distributed uniform random variables on T and are independent of the ij's. 
Noting that the reduced Palm process S i)Q = 0, we obtain the bound by applying (2.22) with 
Pi = P(/j = 1) = T/n and Chen and Xia (2004, Proposition 4.5). ■ 



3 Superposition of thinned dependent point processes 

Assume that q is a measurable retention function on T and X is a point process on T. 
For a realization X(u) of X, we thin its points as follows. For each point of X(u) at a, 
it is retained with probability q(a) and discarded with probability 1 — q(a), independently 
of the other points [see Daley and Vere- Jones (1988), p. 554, for dependent thinning and 
Schuhmacher (2005b) for discussions of more thinning strategies]. The thinned configura- 
tion is denoted as X q {u). For retention functions (ft, • • • ,Qn, let Y17=iX qi be the pro- 
cess arising from the superposition of independent realizations of X qi ,X q2 , . . . ,X qn , that is, 
X' qi ,X' q2 , . . . , X' qn are independent and C (X' q .) = C (X qi ) for i = 1, . . . , n. Fichtner (1977) 
showed that a sequence of such superpositions, obtained from the rows of an infinitesimal ar- 
ray of retention functions, converges to a Poisson process under standard conditions [see also 
Kallenberg (1983, Exercise 8.8)]. Serfozo (1984) presented convergence theorems for sums of 
dependent point processes that are randomly thinned by a two-step procedure which deletes 
each entire point process with a probability and for each retained point process, its points 
are deleted or retained according to another thinning strategy. Necessary and sufficient con- 
ditions are given for a sum of two-step thinned point processes to converge in distribution 
and the limit is shown to be a Cox process [see also Fichtner (1975) and Liese (1980)]. 

For simplicity, we assume that {Sj, i G X} is a locally dependent collection of point 
processes (satisfying [LD1]) on a locally compact metric space V with metric cLq bounded 
by 1. For each point of Sj, we delete the point with probability 1 — p and retain it with 
probability p, independent of the others. The thinned point process is denoted by Sf , i G X, 
and in general, for each point process X, we use X p to stand for its thinned process. Let 
E p = • ^ s before, we define Aj_ as the collection of indices j of the point processes 

Ej which are dependent on Sj, i.e., £j is independent of {Ej, j G Af). 

Theorem 3.1 Let ^ be the mean measure ofEi, ^ = /^j(r) = E(|Hj|) ? i G X and assume 
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that A = Ylieif'i < 00 ■ Then the mean measure ofE p is X p = pY^iei^i an d 

d tv (C(~ p ), Po(X p )) <p(lA±) E^{[|^| + iHiOAi + [\Vi\ + IS,] - 1]|~,|}, (3.1) 



d 2 (£(~ p ),Po(X p )) 
/ 3.5 



2.5 



E 



{[|^| + |Si|]Ai + [|^| + -l]|Si|}, (3.2) 



<WOC(S*),Po(A p )) < pE^UM + |~ t |]A, + + |Si| - 1]|H*|>. 



(3.3) 



Proof. We prove (3.2) only, as the proofs of (3.1) and (3.3) are similar to that of (3.2). By 
conditioning on the configurations, we have, for each Borel set BcT, 

mm] = nmmm = m^u 

which implies that the mean measure of 3? is Af = pfa and hence X p = p J2iex A**- By (2.3) 
and the fact that ^(^1,^2) < |£i| + we obtain 



d 2 (C (E p ),Po(X p )) 





'3.5 


2.5 






+ 1 






2.5 




, pA 




+ 1 




^3.5 


2.5 




1 pA 




+ 1 



-j j(Wl + 1^1 + |3f| + ISfjAf (da) 



{[1^1 + ISfOA^+O^I + lHfl-lJISfl}. 



Since the points are thinned independently, we can condition on the configuration of {Sj, % G 
I}. Noting that for Z ~ Binomial(n,p), < and E[(X - 1)X] = n(n - l)p 2 , we 

obtain 



d 2 (£(3>),Po(A p )) 

3.5 
pX 



4 ex \^ P 



2.5 



This completes the proof of (3.2). ■ 



{[|^| + 13,11a, + + 1341-1113*1}^. 



Remark 3.2 Serfozo (1984, Example 3.6) obtained the rate p for the convergence of a sum 
of thinned point processes to a Poisson process. Theorem 3.1 shows that the rate p is valid 
for all of the three metrics used. 
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4 Superposition of renewal processes 



Viswanathan (1992, page 290) states that if {Hj, 1 < i < n} are independent renewal 
processes on [0, T], each representing the process of calls generated by a subscriber, then 
the total number of calls can be modeled by a Poisson process. In this section, we quantify 
the statement by giving an error bound for Poisson process approximation to the sum of 
independent sparse renewal processes. We begin with a technical lemma. 

Lemma 4.1 Let rj ~ G, & ~ F, i > 1 be independent nonnegative random variables and 
define 

N t = max{n : 77 + £1 H h f„_i < t}, t> 0. 

Then 

G(t) < E(JV t ) < ^^y, (4.1) 
tt~\ / a t \ 2F(rj)E(AQ 2F(t)G(t) , . 

Proof. Let V(t) = E(A r t ). The renewal equation gives 

V(t) = G(t) + y V(i-s)dF(s) 

< G(t) + y(t)F(t), (4.3) 

which implies (4.1). For (4.2), define V 2 (t) = E[A^(A^ + 1)]. Then using the same arguments 
as for proving the renewal equation, 



V 2 (t) = 2V(t) + f V 2 (t - s)dF(s). 
Jo 



This implies V 2 (t) < 2V(t) + V<2(t)-F(t), which in turn implies 

v 2 (t) < -^L. 

Since 

E(7V t 2 ) - E(JVi) = V a (f) - 2V(f) = f V 2 {t - s)dF(s) < V 2 (t)F(t) < ^Jff ' 
(4.2) follows from (4.1). ■ 

Theorem 4.2 Suppose that {Sj, 1 < z < n} are independent renewal processes on [0,T] 
with the first arrival time of having distribution G{ and its inter- arrival time having 
distribution Fi. Let H = J^ILi "» anc ^ ^ ^ e ^ s mean measure. Then 

d (a*) Po(x)) < ^umT)+GmGm 
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Proof. We view a renewal process as a point process whose points occur at the renewal 
times. For a renewal process X with renewal times T\ < r 2 < . . . , we further define X' = 5 T1 . 
Since A, = E(|Sj|), it follows from (2.16) that 

d 2 (£(s),p (A)) < J2 [ x + E y- ... ih 5 .| + 1 ) [A? + E(|s - |2) ~ Xi] 

[A? + E(|S i | a )-A i ]. (4.5) 




* E^l^.| + 1 
However, applying Proposition 4.5 of Chen and Xia (2004) gives 



E — : : < 



and using (4.1), we obtain 

n 

A>^G t (T). 

t=i 

By combining (4.5), (4.1) and (4.2), we obtain (4.4). ■ 

Remark 4.3 If {Sj, 1 < i < n} are independent and identically distributed stationary 
renewal processes on [0, T] with the successive inter-arrival time distribution F, then 

2K K h v ;; - (n-l)(l-F(T)) 2 ' 
where G(t) = f*(l - F(s))ds/ / °°(1 - F(s))ds [see Daley and Vere- Jones (1988), p. 71]. 

Remark 4.4 An application of Theorem 2.1 of Schuhmacher (2005a) to the sum of the 
renewal processes {Hj, 1 < % < n} in Remark 4.3 with the natural partition {{«}, 0, {1, . . . , %— 
1, % + 1, . . . , n}} for each 1 < i < n will give an error bound 

n[F{T) + G(T)\ + 6G(T)(l + ln + n), 

where 9 is a constant. The first term of the bound increases linearly in n and the bound is 
clearly not as sharp as the bound in Remark 4.3. 

Since the thinned process X p of a renewal process X with mean measure fi is still a 
renewal process [see Daley and Vere- Jones (1988, pp. 75-76)] with mean measure /x p = p/j, 
[see the proof of Theorem 3.1], a repetition of the proof of Theorem 4.2 yields the following 
proposition. 

Proposition 4.5 Suppose that {Hj, 1 < i < n} are independent renewal processes on [0, T\ 
with the first arrival time of Sj having distribution Gi and its inter- arrival time having 
distribution Fj. Let Hf be the thinned point process obtained from Ej by deleting each point 
with probability 1 — p and retaining it with probability p, independently of the other points. 
Let E p = Yli=i an d ^ P be its mean measure. Then 



C / 2 (£(E^),Po(A^))< 



(ELi Gi(T) - max, G 3 (T)) (1 - F z (T)y 



13 



Acknowledgements 

This research was partially supported by Grant C-389-000-010-101 at the National Uni- 
versity of Singapore (LC), the ARC Centre of Excellence for Mathematics and Statistics of 
Complex Systems (AX) and the Belz Fund of the University of Melbourne (LC). 

References 

[1] Banis, R. (1975). The convergence of sums of dependent point processes to Poisson 
processes. (Russian. Lithuanian, English summary) Litovsk. Mat. Sb. 15, 11-23, 
223. 

[2] Banis, R. (1985). A Poisson limit theorem for rare events of a discrete random 
field. (Russian. English, Lithuanian summary) Litovsk. Mat. Sb. 25, 3-8. 

[3] Banys, R. (1980). On superpositions of random measures and point processes. 
Mathematical statistics and probability theory (Proc. Sixth Internat. Conf., Wisla, 
1978), 26-37, Lecture Notes in Statist., 2, Springer, New York-Berlin. 

[4] Barbour, A. D. (1988). Stein's method and Poisson process convergence. J. Appl. 
Probab. 25 (A), 175-184. 

[5] Barbour, A. D. and Brown, T. C. (1992). Stein's method and point process approx- 
imation. Stochastic Processes and their Applications 43, 9-31. 

[6] Barbour, A. D. and Eagleson, G. K. (1983). Poisson approximation for some statis- 
tics based on exchangeable trials. Adv. Appl. Prob. 15, 585-600. 

[7] Barbour, A. D., Hoist, L. and Janson, S. (1992). Poisson Approximation. Oxford 
Univ. Press. 

[8] Barbour, A. D. and Xia, A. (2006). Normal approximation for random sums. Adv. 
Appl. Prob. 38, 693-728. 

[9] Brown, T. C. (1978). A martingale approach to the Poisson convergence of simple 
point processes. Ann. Probab. 6, 615-628. 

[10] Brown, T. C. (1979). Position dependent and stochastic thinning of point processes. 
Stochastic Processes and their Applications 9, 189-193. 

[11] Brown, T. C, Weinberg, G. V. and Xia, A. (2000). Removing logarithms from 
Poisson process error bounds. Stochastic Processes and their Applications 87, 149- 
165. 

[12] Brown, T. C. and Xia, A. (1995). On Metrics in Point Process Approximation. 
Stochastics and Stochastics Reports 52, 247-263. 

[13] Chen, L. H. Y. and Shao, Q. M. (2004). Normal approximation under local depen- 
dence. Ann. Probab. 32, 1985-2028. 



14 



[14] Chen, L. H. Y. and Xia, A. (2004). Stein's method. Palm theory and Poisson 
process approximation. Ann. Probab. 32, 2545-2569. 

[15] Cinlar, E. (1972). Superposition of point processes. Stochastic point processes: 
statistical analysis, theory, and applications (Conf., IBM Res. Center, Yorktown 
Heights, NY., 1971), 549-606. Wiley-Interscience, New York. 

[16] Daley, D. J. and Vere- Jones, D. (1988). An Introduction to the Theory of Point 
Processes. Springer- Verlag. 

[17] Fichtner, K. (1975). Schwache Konvergenz von unabhdngigen Uberlagerungen 
verdunnte zufdlliger Punkfolgen. Math. Nachr. 66, 333-341. 

[18] Fichtner, K. (1977). Poissonsche zufdllige Punktfolgen und ortsabhdngige 
Verdiinnungen. (German) Transactions of the Seventh Prague Conference on In- 
formation Theory, Statistical Decision Functions, Random Processes and of the 
Eighth European Meeting of Statisticians (Tech. Univ. Prague, Prague, 1974), Vol. 
A, pp. 123-133. 

[19] Goldman, J. R. (1967). Stochastic point processes: limit theorems. Ann. Math. 
Statist. 38, 771-779. 

[20] Goldstein, L. and Rinott, Y. (1996). Multivariate normal approximations by Stein's 
method and size bias couplings. J. Appl. Probab. 33, 1-17. 

[21] Goldstein, L. and Xia, A. (2006). Zero Biasing and a Discrete Central Limit The- 
orem. Ann. Probab. 34, 1782-1806. 

[22] Grigelionis, B. (1963). On the convergence of sums of random step processes to a 
Poisson process. Theory Probab. Appl. 8, 177-182. 

[23] Jagers, P. (1972). On the weak convergence of superpositions of point processes. Z. 
Wahrsch. Verw. Gebiete 22, 1-7. 

[24] Kallenberg, O. (1975). Limits of compound and thinned point processes. J. Appl. 
Probab. 12, 269-278. 

[25] Kallenberg, O. (1983). Random Measures. Academic Press, London. 

[26] Liese, F. (1980). Uberlagerung verdiinnter und schwach abhdngiger Punktprozesse. 
(German) Math. Nachr. 95, 177-186. 

[27] Rachev, S. T. (1991). Probability metrics and the Stability of Stochastic Models. 
John Wiley & Sons . 

[28] Serfozo, R. (1984). Thinning of cluster processes: convergence of sums of thinned 
point processes. Math. Oper. Res. 9, 522-533. 

[29] Schuhmacher, D. (2005a). Distance estimates for dependent superpositions of point 
processes. Stochastic Processes and their Applications 115, 1819-1837. 



15 



[30] Schuhmacher, D. (2005b). Distance estimates for Poisson process approximations 
of dependent thinnings. Electron. J. Probab. 10, 165-201 (electronic). 

[31] Viswanathan, T. (1992). Telecommunication switching systems and networks. 
Prentice-Hall of India. 

[32] Xia, A. (2005). Stein's method and Poisson process approximation. In: An Intro- 
duction to Stein's Method, Eds. A. D. Barbour and L. H. Y. Chen, World Scientific 
Press, Singapore, 115-181. 



16 



